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For short periods of time during June and July of 1963 the U.S. Navy transmitters 
located at Jim Creek, Wash. (NPG), and Balboa, Panama (NBA), simultaneously trans- 
mitted a frequency stabilized signal of 18 kc/s. The phase and relative amplitude of these 
signals were monitored at Boulder, Colo.; College, Alaska; Maui, Hawaii; and Tucuman, 
Argentina. A semiempirical method by which the mean relative phase velocity can be 
calculated from the measurements of phase made at these four recording sites is demon- 
strated. The values of the mean relative phase velocity at 18 kc/s which are obtained for 
daytime and nighttime propagation conditions are respectively, (V v /c) d = 0.998 and 
(V p /c) n = 0.995. 



1. Introduction 

Early experimental work [Al'pert and Borodina, 
1959; Jean, Taylor, and Wait, 1960] on the velocity 
of propagation of VLF radio waves utilizing the 
electromagnetic energy radiated from lightning dis- 
charges seemed to indicate for frequencies lower than 
about 20 kc/s that the phase velocity in the earth- 
ionosphere waveguide would be appreciably higher 
than the velocity of light, c in vacuum. Taken 
together the above cited work gives essentially the 
same result for both day [Al'pert and Borodina, 1959] 
and night [Jean et al., 1960] propagation conditions. 
However, more recent VLF propagation measure- 
ments, obtained by means of the Radux-Omega 
navigation system [Casselman et al., 1959], have 
shown [Wait, 1961] that during nighttime propaga- 
tion conditions the phase velocity is apparently 
lower than c for frequencies greater than about 9 
kc/s, and under daytime conditions lower than c for 
frequencies greater than about 14 kc/s. 

Knowledge of the phase velocity in the earth- 
ionosphere waveguide is an important requirement 
for the proper interpretation of the phase variations 
observed in VLF propagation data and the observed 
systematic errors in phase-comparison types of 
VLF navigation aids. Thus, it is the purpose of this 
paper to present a method by which an estimate of 
the phase velocity can be obtained and, as an appli- 
cation of this method, to give the results obtained at 
a frequency of 18 kc/s. 

2. Method of Measurement and Basic 
Theory 

For short periods of time during June and July of 
1963, the U.S. Navy transmitters NBA at Balboa, 
Panama, and NPG at Jim Creek, Wash., simul- 
taneously transmitted phase stabilized signals at a 
frequency of 18 kc/s [F. M. Malone, private communi- 
cation, 1963]. Measurements of the phase of the 



two signals were made using the transmitter- 
receiver network illustrated in figure 1. The trans- 
mitter-receiver pairs required to make the measure- 
ments, from which a calculation of relative phase 
velocity can be obtained, is illustrated diagrammat- 
ically in figure 2. The symbols A\ and R 2 represent 
the receivers, while 2\ represents the transmitter 
NPG at Jim Creek, Wash., and T 2 the transmitter 
NBA at Balboa, Panama. The propagation path 
lengths are D if D 2y D 3 , and D A . The received phases 
of NBA and NPG at R x (Boulder) and R 2 (Maui) are 
respectively ($i, $ 2 ) and (<£ 3 , $ 4 ). It was possible to 
measure the received phase of each signal at four 
monitoring sites: College, Alaska; Boulder, Colo.; 
Maui, Hawaii; and Tucuman, Argentina. These 
four receiving sites were combined to give six receiver 
pairs recording the NBA and NPG transmissions. 




Figure 1. Geographic location of the VLF transmitter-receiver 
networks showing the propagation paths used in this study of 
relative phase velocity. 

Each dash represents approximately 10 wavelengths of the 18 kc/s transmitter 
frequency. 
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Figure 2. Diagrammatic representation of an individual trans- 
mitter-receiver network, upon which the calculation of a 
relative phase velocity is based. 

The vector diagrams illustrate the cycle ambiguity arising at each receiving 
site in the calculation of the measured phase difference. 



The data pairs used for this analysis are: Maui- 
Tucuman, College-Boulder, College-Maui, College- 
Tucuman, Boulder-Maui, and Boulder-Tucuman. 
It will be shown in this paper that by assuming a 
uniform, isotropic, homogeneous ionosphere for the 
portion of the waveguide containing the two trans- 
mitters and four receiving sites, an effective value of 
the relative phase velocity, V p /c } can be estimated 
from the data obtained from these pairs of recording- 
sites. 

; A typical record, obtained at Boulder, Colo., is 
shown in figure 3. In the time period from approx- 
imately 1800 UT to 1900 IJT, both NBA and NPG 
transmitted on 18 kc/s. The trace at the top of the 
record is the relative amplitude; the arrow on the 
right gives the direction of an amplitude increase. 
The trace at the bottom of the record is the phase; the 




HOURS, UT 

Figure 3. Copy of an original record of the phase and ampli- 
tude of the signals NBA (18 kc/s, Balboa, Panama) and 
NPG (18 kc/s, Jim Creek, Wash.) as recorded at Boulder, 
Colo. 

arrow at the left gives the direction of a phase 
advance. The phase of NBA alone is indicated in 
the time period before 1800 UT. The phase of NPG 
alone occurs in the time period from 1807 UT to 1810 
UT when NBA is regularly not transmitting. Thus 
in these time periods at Boulder on 11 July 1963, the 
received phase ($i) of NBA is 234° and the received 
phase ($ 2 ) of NPG is 327°. For purpose of illus- 
tration these values are shown in the vector diagram 
of figure 2. Similar records were taken at the other 
three recording sites, and the measured phase values 
of NPG and NBA from all of the usable records are 
tabulated (table 1). In table 1, the location of each 
recording site is given together with the date and 
measurement period in Universal Time. The data 
were separated according to solar controlled propa- 
gation conditions; i.e., either the propagation path 
was totally dark (all dark) or totally sunlit (all light) . 
Each record gave a value of phase for NPG and NBA 
with respect to the receiver reference oscillator. Due 
to the electronic characteristics of the equipment, a 
possible cycle ambiguity (fig. 2) has to be taken into 
account when the phase difference between these two 
values is calculated. 







Table 1. 


Measured VLF phase data 










Date 
1963 


Measure- 
ment, 
period 
UT 


Phase 


Propagation 


Location of recording site 


condition 


Boulder 


College 


Maui 


Tucuman 




NPG 


NBA 


NPG 


NBA 


NPG 


NBA 


NPG 


NBA 


Path all dark 


27 June 

28 June 

29 June 

30 June 
5 July 

26 June 
26 June 
9 July 

11 July 

12 July 
14 July 


0600-0800 
0600-0800 
0600-0800 
0600-0800 
0600-0800 

1300-1500 
2100-2400 
1800-2100 
1800-2100 
1800-2100 
2100-2400 


deg 
167 
274 
216 
162 
176 

270 
263 


deg 
158 
291 
338 
320 
198 

342 
346 


deg 
248 
241 
198 


deg 
155 
173 
234 


deg 
95 
36 
324 


deg 
36 
353 
36 


deg 
223 
295 
192 
112 
263 


deg 
97 
220 
135 
119 












144 


Path all light 


288 
288 


65 
76 










238 

207 
274 
155 
111 


317 
112 
144 
162 

5 


329 
173 
252 

25 


202 
230 




327 

194 


234 

277 






280 








220 








245 
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In the analysis of the data we make use of the 
expression [Wait, 1961] that relates the phase 
velocity (V p ) to the total phase path between the 
transmitter and receiver in the earth-ionosphere 
waveguide : 



Vr, 



= 2ttM+0, 



where co is the angular wave frequency, D is the path 
distance, M is the integral number of wavelengths 
in the propagation path, and is the additional 
fractional part of a cycle in the phase path. Thus 
we have for the propagation path system shown in 
figure 2 the four equations: 






=2wM i +ct> i , 



i=l,2, 3,4. 



Combining these four equations, and noting that 
the measured phase differences are given by 
$i— $2=01— </>2+Y and $ 3 — $4 = 03 — #4+7> results 
in two expressions: one for the phase difference 
measured at receiver R h 

($ 1 _$ 2 ) = _2 7 r(M 1 -M 2 ) + (co/y,)(A-A)+7 ? (1) 
and one for the phase difference measured at receiver 

xl2> 

(*4-* 4 )= -2tt(M 3 -M 4 ) + (co/F p )(A-A)+7, (2) 

where y is the phase offset between the two trans- 
mitters. Equations (1) and (2) are then combined 
and solved for the phase velocity: 



V^MtiDt) 



o[(A-A)-(A-A)] 



($i-$ 2 )-($3-$4)+27r[(M 1 -M 2 )-(M 3 -M 4 )] 
By letting V p =c we obtain an auxiliary equation: 

c(0;,m;,a) 

= co[(A-A)-(A-A)] 

(0;-02)-(03-0O+27t[(m;-mo-(M^-mi)]' 

where M' t and 0- are respectively the calculated 
integral number of wavelengths, and additional 
fractional part of a cycle in the propagation path. 
These primed quantities M\ and 0- are calculated 
from the four equations: 

^=2ttM;-+0;-, { = 1,2,3,4, 

where the geographic distances D L are computed, 
using the best available geographic coordinates, by 
the method of Lambert [1942] in which the earth is 
treated as an ellipsoid of revolution. Since Lam- 
bert's [1942] method gives a value of path distance 
which appears to differ from the true distance by 
less than one part in 200,000 [Sitterly and Pierce, 



1944] even when the distance is nearly half a cir- 
cumference, these values of D t are assumed to be 
within one wavelength (at 18 kc/s, approximately 
1G km) of the true path length. 

Finally, by taking the ratio V p /c we obtain the 
relation for the relative phase velocity in terms of 
measured and calculated phase path parameters: 

F / ,_(0;-0O-(0^-0;)+2 7 r|(.u;-^)-(A^-Mi)] 



C ($1 — $ 2 ) — ($3 



-$ 4 )+2tt[(M 1 -M 2 )-(M 3 -M 4 )] 

(3) 



Using (3), the calculation of V p /c can be carried out 
by taking, as a starting point, the values of M t 
equal to the calculated M-; then the M.- t are in- 
dividually varied by integral numbers of wavelengths. 
Thus by a method of successive approximations a 
set of independent values is obtained for each 
receiver pair of measurements. Taken together 
these independent sets of values uniquely determine 
a value of V p /c which is common to each set. Then 
the average of these common values of relative 
phase velocity can be taken as the best estimate. 

3. Discussion and Analysis of the Data 

The eight VLF propagation paths, determined by 
the two transmitters and four receiving sites, give a 
total of six path systems for which calculations of 
the relative phase velocity V p /c can be carried out 
using (3). The cycle ambiguity illustrated in figure 2 
and the variation of the M t by integral numbers 
result in four independent trial values of the ratio 
Vp/c for each receiver pair. If any one of the four 
values is significant, this value should repeat itself 
at different times and in each set of calculations for 
(Mich of the receiver pairs (Maui-Tucuman, College- 
Boulder, College-Maui, College-Tucuman, Boulder- 
Maui, and Boulder-Tucuman). There were 11 
periods consisting of 6 days and 5 nights when at 
least 2 of the 4 receiving sites obtained usable data. 
The data obtained during these 1 1 recording periods 
were analyzed in the manner outlined above and 
yielded a total of 140 calculated values of V p /c. 
Sixty of these values were obtained for daytime 
propagation conditions, and 80 were obtained for 
nighttime conditions. The results of this analysis 
for nighttime conditions are shown in figure 4. In 
this figure all of the calculated trial values of relative 
phase velocity for each measurement during night- 
time propagation conditions are grouped according 
to receiver pairs as indicated. Thus each receiver 
pair shows columns of four possible values of the 
relative phase velocity where each column represents 
an independent experimental measurement for that 
receiver pair. There were four receiving sites which 
were combined to give six possible receiving pairs. 
That this combination constitutes an inherent check 
on the internal consistency of the data is made 
evident by the fact that in each set of 4 calculations 
(each column) only one value is found to be common 
to all 20 sets (columns). Although a considerable 
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RECEIVER PAIRS 



Figure 4. Calculated trial values of relative phase velocity for 
nighttime propagation conditions grouped according to receiver 
pairs. 

Each column of four values represents an individual experimental measure- 
ment. The most probable, value of (V p /c) n is that which is common to each 
group and occurs between the dashed lines. 

spread in each set of four trial values is evident, a 
value near 0.995 consistently appears in each 
grouping. These 20 values are set apart by dotted 
lines in figure 4. The mean of these common values 
is 0.995. For daytime conditions the analysis indi- 
cates a value near 0.998 as the common value of 
relative phase velocity (fig. 5). The mean of these 
daytime values is 0.998. 

A standard deviation of the mean for the above 
relative phase velocities was not given because (a) 
some of the data for the daytime case was rejected, 
and (b) only two in the six possible receiver pairs 
can be considered as being strictly independent. 
However, for the 20 values of the nighttime case 
the sample standard deviation is 0.0003, giving 
0.0001 as an estimate for the standard deviation of 
the mean even if the effective number of strictly 
independent values (of the 20) is as small as 5. In 
the daytime case, including the rejected 4 points to 
give a total of 15 values, the sample standard 
deviation is 0.0023, giving 0.0009 as an estimate for 
the standard deviation of the mean where the effec- 
tive number of strictly independent values (of the 
15) is as small as 7. However, if in the daytime 
case we use only the 11 accepted values lying between 
the dotted lines of figure 5, the sample standard 
deviation is 0.0004, giving 0.0002 as an estimate for 
the standard deviation of the mean even if the 
effective number of strictly independent values (of 
the 11) is as small as 7. In this last case it is clear 
that this is a lower limit since it is possible that the 
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Figure 5. Calculated trial values of relative phase velocity for 
daytime propagation conditions grouped according to receiver 
pairs. 

Each column of four values represents an individual experimental measure- 
ment. The most probable value of (V p /c)d is that which is common to each 
group and occurs between the dashed lines. 

random errors, as evidenced by the four rejected 
values, may have exceeded the inherent ambiguities. 

In the 20 sets of calculations for the nighttime rela- 
tive phase velocity, a common value occurred 20 
times. However, in the 15 sets of trial values for 
the daytime relative phase velocity, (V p /c) d , a com- 
mon value occurred only 11 times. The fact that 
the common value did not occur as consistently in 
the daytime calculations as in the nighttime, is diffi- 
cult to explain. The four values in the daytime cal- 
culations which did not fall into the general pattern 
arose from the pairing of the low magnetic latitude 
data from Maui, Hawaii with the higher magnetic 
latitude data from Boulder, Colo., and College, 
Alaska. This might point to a possible latitude vari- 
ation in the height of the daytime ionosphere [Chil- 
ton et al., 1964]. Nevertheless, other parameters 
such as the direction of propagation may play an 
important role and should therefore be considered. 
However, the limited amount of data presently avail- 
able is not sufficient for the resolution of this problem. 

The above values of experimentally deduced rela- 
tive phase velocities are in good agreement with the 
theory of the first order mode for VLF propagation 
in a curved earth-ionosphere waveguide [Wait, 1963]. 
Specifically, theoretical calculations by Spies and 
Wait [1961] show that for a reasonable model of the 
ionosphere and for a perfectly conducting ground the 
relative phase velocities obtained above correspond 
to effective heights of reflection of 70 km during the 
day, and 95 km at night. Since the propagation 
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paths used in this study cover a large variation in 
latitude and longitude these values of effective height 
should only be taken as an approximation. 'Phis is 
emphasized by the fact that the difference in calcu- 
lated diurnal height change between the NBA- 
Boulder path and the NBA-Tacuman path is on the 
order of 5 km [Chilton et al., 1964] even though the 
two paths are essentially the same length. 

It is difficult to compare the above phase velocities 
with the earlier broadband experimental results 
[Jean et al., 1960; APpert and Borodina, 1959] due 
to the fact, not specifically stated in either paper, 
that in their analysis the phase velocities at higher 
frequencies were assumed and taken to be equal to 
the velocity of light. On the other hand, our values 
at 18 kc/s compare favorably with extrapolated 
values from the later work [Wait, 1961] done in the 
frequency range from 8 kc/s to 16 kc/s. 

4. Conclusions 

This paper shows that for a frequency of 18 kc/s 
and for that portion of the earth-ionosphere wave- 
guide covered by the transmitter-receiver network 
shown in figure 1, the mean effective relative phase 
velocities are (V p /c) d =0.998 by day and (V p /c) n = 
0.995 by night. 

Although the experimental technique and method 
of analysis of the data are relatively uncomplicated, 
they yield results which are strongly dependent on 
having enough data to minimize the ambiguities and 
experimental errors which are inherent in the meas- 
urements. There is some evidence in the observa- 
tions reported here which might suggest that the 
daytime phase velocity depends on path latitude and 
possibly direction of propagation. However, fur- 
ther observations on more paths and at different fre- 
quencies are required to verify this conclusion. 
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Several comments are in order on the phase veloci- 
ties derived by Steele and Chilton [1964]. 

1. The NPG and NBA to College paths were not 
totally dark, as claimed, between 0600 and 0800 UT 
in late June and early July. In fact, College was 
actually in sunlight at those times. By the end of 
their observation period on July 13, only roughly 25 
percent of the D region at 0600 was dark on the NPG- 
College path, and 75 percent on the NBA-College path; 
by 0800 the amounts that were dark increased to about 
60 and 90 percent, respectively. The inclusion of the 
College data in their computation of the nighttime 
phase velocities therefore seems questionable. 

2. Nighttime second mode propagation effects can- 
not be ignored on any of the paths discussed, partic- 
ularly the NPG to Boulder and College paths. Both 
of these paths are so short that the second mode may 
well predominate at night [Watt and Croghan, 1964]. 
In addition, the observed nighttime phase and ampli- 
tude effects observed on various longer VLF paths, 
including some of the paths Steele and Chilton discuss, 
can only be consistently explained by a mode-inter- 
ference model in which the first and second modes are 
of roughly comparable amplitude on 4000 to 5000 km 
paths [Bates and Albee, 1965]. Crombie's [1964] 
apparently correct explanation of the sunrise fading 
effect relies upon a strong nighttime second mode 



propagating to distances well over 7500 km; quantita- 
tive bounds on the second mode amplitude can be ob- 
tained from his model. 

The equations used by Steele and Chilton are valid 
only for single-mode propagation. Therefore, the 
nighttime phase velocities they found should be recom- 
puted with phase interference effects included. 

3. It is not specifically stated whether the most 
probable values of phase velocity they give corre- 
spond to the same day and night values of the M\ for 
each path. This would appear to be a test for con- 
sistency. Also, the values of Mi found for each path 
would seem to be an important parameter of VLF 
propagation. 
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In the following we reply to Bates' comments in the 
order of their presentation. 

1. Our calculations show that on 28 June 1963 the 
average solar zenith angle (x) over the NPG-College 
path varies from 90° to 94° during the measurement 
period (0600 to 0800 UT), and for the NBA-College 
path varies from 90° to 94° during the measurement 
is from 94° to 96°. On the basis of these calculations 
we have called the path totally dark when the average 
zenith angle is greater than 90°. However, we agree 
with Bates that portions of the two College paths were 
weakly illuminated, but the paths can be considered 
to be more nearly representative of nighttime con- 
ditions than of daytime conditions. Nevertheless, 
the inclusion of the College data in the computations 
of the nighttime phase velocities might seem question- 
able. Let us assume for the moment that the College, 
Alaska, data are marginal. Eliminating these data 
from the analysis does not change the result, i.e., the 
calculated average value of the phase velocity remains 
the same. 

2. Whether or not the nighttime second mode 
propagation effects can be ignored depends on the 
relative attenuation (a) of the second-order mode 
compared to the first-order mode, and the magnitude 
of the excitation factor (A) for the second-order mode 
compared to the first-order mode. For frequencies 
greater than 20 kc/s both theory [Wait and Spies, 
1964] and experiment [Watt and Croghan, 1964] 
indicate that there may very well be an appreciable 
second-order mode effect. However, for frequencies 
less than 20 kc/s both theory [Wait and Spies, 1964] 
and experiment [Wait, 1961; Watt and Croghan, 1964] 
indicate that the effect of the second-order mode is 
small if not negligible at distances greater than 3000 
km. In point of fact, our measurements were made 
at 18 kc/s, and all but two of the eight propagation 
paths had distances greater than 4000 km. The 
results obtained using the two shorter paths do not 
appear to be inconsistent with those obtained from the 
rest of the data. In regard to Bates' interpretation 
of Crombie's [1964] explanation of the sunrise fading 
effects, this theory does not necessarily require a 
"strong nighttime second mode propagating to dis- 



tances well over 7500 km," but the observed results 
could possibly be due to additive effects of mode con- 
version in a distributed sunrise region [Crombie, 
1965 private communications]. In addition, it should 
be remembered that Crombie's published results 
refer to paths where propagation was to the east, 
whereas the paths in the paper under discussion were 
mainly for propagation to the west, in which case 
attenuation of the second-order mode is considerably 
increased relative to that of the first-order mode. 

3. The Mi are the integral number of wavelengths 
in the respective propagation paths. In the equation 
(3), used to calculate the relative phase velocity, they 
are additively summed: [(Mi - M 2 ) — (M 3 — M 4 )]. We 
did not attempt to evaluate an absolute Mi for each 
path, because (a) the Mi are not initially known or meas- 
ured but are obtained as a first approximation by 
calculation assuming the velocity of light as the phase •> 
velocity, and (b) there is a cycle ambiguity, inherent in 
each measuring equipment (see fig. 2), which could 
amount to a wavelength. It was therefore necessary 
to vary the estimated M t , and thus the sum, by integral 
numbers of wavelengths. However, we agree with 
Bates that the Mi are important parameters in VLF 
propagation and, if one assumes our values of phase 
velocity at 18 kc/s, they can be calculated for each 
path by evaluating the expression [Wait, 1961] that 
relates the phase velocity to the total phase path be- 
tween the transmitter and receiver in the earth-iono- 
sphere waveguide. 
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